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^sO . We investigate the magnetocrystalline anisotropy (MCA) energy of tetragonal distorted FeCo alloys depend- 

ing on the degree of order by first-principles electronic structure calculation combined with the coherent 
potential approximation. The obtained results indicate that the MCA energy of FeCo alloys strongly depends 
on the degree of order under optimal conditions, where the axial ratio of the bet structure is 1.25 and the 
composition is Feo.sCoo.s- We find that the modification of the electronic structure resulting from electron 
scattering by chemical disorder has a considerable influence on the MCA under these conditions. 
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In recent elemental strategies, it has been essential for hard magnetic materials 
applied to magnetic devices and permanent magnets to be free from rare-metal and rare- 
earth elements. Given this situation, Burkert et al. proposed that the FeCo alloy can be 



^ ■ expected to exhibit a giant magnetocrystalline anisotropy (MCA) energy under certain 

OO . conditions by the first-principles calculation using the virtual crystal approximation 

1 (VCA) to treat the binary alloy system. 1 In their paper, tetragonal distorted Fe\- x Co x 

alloys exhibit a large uniaxial magnetic anisotropy constant K u , as high as about 700 
- 800 /ieV/atom, which is sufficiently comparable to the K n value of a series of Lln- 
type alloys such as FePt, CoPt, and FePd, 2 under the conditions that the axial ratio 
^ ■ between the a and c axes, c/a, is 1.20-1.25 and that the composition x is 0.5-0.6. This 

b ■ 

result has stimulated considerable expectation for their use as novel materials with 
high coercivity; thus, many experimental trials have been performed to realize such a 
property. 3-7 Under the optimal conditions, the magnetic easy axis coincides with the 
c-direction in agreement with the theoretical prediction; however, satisfactory results 
have not yet been obtained at a quantitative level. 

Recently, Neise et al. 8 have performed a first-principles investigation on the MCA 
energy of FeCo alloys with a chemical disordered structure by employing both the VCA 
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and the supercell method. When using the VCA, they obtained results similar to those 
of Burkert et al.; 1 however, for the supercell structure reflecting chemical disorder, 
they obtained values of the MCA energy smaller by factors of 1.5-3.0. Therefore, we 
can consider that one of the reasons for the disagreement with the experiments may 
be the chemical disorder. Since the VCA is generally the lowest order approximation 
in terms of a random configuration of atoms, it may be inappropriate to calculate 
the MCA energy in FeCo owing to the insufficient treatment of the chemical disorder. 9 
However, the physical explanation for this discrepancy and the influence of the chemical 
disorder on the MCA in FeCo are not yet understood. In this work, we investigate the 
MCA in disordered FeCo alloys by means of the coherent potential approximation 
(CPA), 10,11 which provides us with a physical picture of the effects of chemical disorder 
in a transparent manner. Employing the CPA has the advantage of reflecting electron 
scattering by the chemical disorder, which has not been considered in the VCA and the 
supercell approach. In addition, it allows us to continuously change the order parameter 
S of the crystal, 12-14 a parameter that is measurable in experiments. 

Electronic structure calculations are performed by using the tight-binding linear 
muffin-tin orbital (TB-LMTO) method under the local spin-density functional approx- 
imation. 15-17 Chemical disorder is taken into account by means of both the CPA and 
the lowest order VCA to compare the obtained results. In the former approximation, 
the averaged potential function (P), which is used in the TB-LMTO method and re- 
flects the random distribution of Fe and Co atoms, is determined via the iterative CPA 
condition for single-site t-matrices. In the latter, (P) is determined only by the sim- 
ple weighted average of the potential functions of both atoms in accordance with the 
composition ratio, 16 so there is no electron scattering. In both methods, we take advan- 
tage of the Green's function technique, which needs an infinitesimal imaginary value 
reflecting causality; we set this value to 2.0 mRy in this work. The MCA energies are 
calculated based on the force theorem with inclusion of the spin-orbit interaction after 
a self-consistent calculation. The lattice constants are determined by the volumes of 
the unit cell of the experimental data of the bulk Fe-Co alloys 18 and the optional axial 
ratio c/a. The number of /c-points used for MCA energy calculation is about 50 3 in the 
full Brillouin zone. 

First, we examine the K n value of completely disordered Fe .5Co .5 (S = 0.0) as a 
function of the axial ratio c/a in Fig. 1. In both the CPA and the VCA, we can measure 
the peak of K u around c/a = 1.25. We note that K u is zero for c/a = 1.0 and \/2 
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owing to the cubic symmetry. This behavior is qualitatively consistent with the results 
of Burkert et al., 1 although the value of K n calculated using the VCA is almost half 
at the peak position. This discrepancy may be attributed to the difference in the band 
calculation method; Burkert et al. 1 used the full-potential LMTO method, which can 
provide more accurate results than the present method if one calculates the electronic 
structure of a perfectly ordered crystal, i.e., S = 1.0. Here, it should be emphasized 
that the value of K n calculated using the CPA is further reduced compared with that 
calculated using the VCA, especially around c/a = 1.25. This behavior is consistent 
with the results obtained by Neise et al. using the supercell method. 8 

Next, let us look at the degree of order dependence of K u for c/a = 1.25 in Feo.sCoo.5 
in Fig. 2. Here, we note that there are two equivalent sites in the primitive cell of the 
tetragonal lattice. We set the composition at the body center site and the corner one in 
the primitive cell to Fei-^Co,, and Fe^Coi^ (0 < rj < 0.5), respectively; the long-range 
order parameter of the crystal is defined by S = 1 — 2t] (0.0 < S < 1.0). 14 For S = 1.0, 
the lattice is a perfectly ordered structure; therefore, the results obtained using the 
VCA and CPA are invariant. In Fig. 2, there is a large difference between these results 
for S < 0.8; therefore, we can consider that the values obtained using the lowest order 
VCA possibly overestimate the uniaxial anisotropy constant of FeCo compared with the 
values obtained using the CPA, which is a more proper way of dealing with a disordered 
alloy system. 11 Therefore, this result reveals that the giant MCA expected from FeCo 
is not present even if there is a slight disorder in the crystal. 

To gain insight into the physical feature giving rise to the different behavior of 
the degree of order dependence of the MCA calculated using the VCA and CPA, we 
focus on the partial Bloch spectral functions of the d x 2_ y 2 and d xy orbital components 
in the minority-spin states at the T-point plotted in Fig. 3. According to Burkert et 
al., these two orbital components have a major contribution to the MCA energy of 
tetragonal distorted FeCo alloys such that the magnetic easy axis coincides with the 
c-direction. This can be understood by expressing the energy variation arising within 
the second-order perturbation in terms of the spin-orbit interaction "Hso/ 9 ' 20 



where Ekn\. denotes the eigenvalue of the nonperturbative eigenstate \kn\) in the 
minority-spin state. Thus, the projected components {d xy \kni) ■ \d xy ) and {d x 2_ y 2 \kni) ■ 
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\d x 2_ y 2) mainly contribute to SE and MCA energy. The restriction of the summation 
is applied to the occupied (unoccupied) state below (above) the Fermi level; therefore, 
eq. (1) indicates that the MCA energy becomes large when the Fermi level lies be- 
tween d x 2_ y 2 and d xy states. As shown in Fig. 3, the partial Bloch spectral function 
of FeCo (x = 0.5, c/a = 1.25) for (a) ordered structure (S = 1.0) and (b) disordered 
one (S = 0.0) calculated using the VCA certainly supports this picture. The peaks of 
the spectra of the occupied d x 2_ y 2 and unoccupied d xy states are coupled to each other 
through the spin-orbit interaction, thereby inducing a large MCA. The MCA energy is 
more enhanced for S = 1.0, since the energy separation between the occupied d x 2_ y 2 
and unoccupied d xy states is smaller than that for S = 0.0 calculated using the VCA. 
This is one of the reasons why the uniaxial anisotropy constants increase with S in Fig. 
2. 

However, the spectra for S = 0.0 obtained via the CPA shown in Fig. 3(c) also seem 
to meet this condition, although it should be noted here that in Fig. 3(c), the shapes 
of these spectra around the Fermi level are modified owing to electron scattering by 
the random arrangement of atoms. The width of the spectra increases compared with 
that obtained using the VCA as shown in Fig. 3(b). Some portions of the spectrum 
lie in the occupied state and others lie in the unoccupied states because of the wide 
energy distribution of each spectrum around the Fermi level. When we evaluate the 
energy variation 5E that results from mixing d xy and d x 2_ y 2 states in the minority-spin 
state through the spin-orbit interaction by using the nonperturbative Green's function 
obtained using the CPA and VCA based on the approach of Solovyev et al., 21 we 
confirm that 5E calculated using the CPA is about one-half of that found using the 
VCA; 22 therefore, this modification of the spectra contributes to decreasing the MCA 
energy in the bet FeCo alloy in Figs. 1 and 2. Since in the CPA calculation, the energy 
widths of these two states increase with decreasing S, the K u values from the CPA 
decrease much faster than those from the VCA when S decreases. Therefore, the giant 
MCA energy of the bet FeCo alloy can be realized under the quite strict condition 
concerning the d xy and d x 2_ y 2 states at the T-point. 

In summary, we investigated the dependence of MCA energy on the degree of order 
of FeCo alloys by means of a first-principles calculation combined with the CPA. The 
results indicate that the MCA energy of FeCo alloys strongly depends on the degree 
of order under the optimal condition of the axial ratio c/a = 1.25. This is because the 
modification of the spectral function resulting from electron scattering in the chemical 
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disordered structure has a considerable influence on the MCA energy under this condi- 
tion, and so the K n value of the FeCo alloy obtained using the CPA is about one-half 
of that obtained using the VCA; therefore, electron scattering has a significant influ- 
ence on the MCA energy, which gives rise to the large difference between the VCA and 
CPA calculations in the small-S" region. Finally, we should stress that the high degree 
of order, in addition to the large distortion of the lattice, is a necessary condition for 
obtaining the large MCA from FeCo in experiments. 
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Fe 5 Co 5 (S=0.0) 




Fig. 1. Uniaxial anisotropy constant K u of FeCo for S = 0.0 as a function of the axial ratio c/a. 
The results obtained from the CPA and VCA are denoted by the open circles and filled squares, 
respectively. 
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Fig. 2. Uniaxial anisotropy constant K u as a function of degree of order S of FeCo with c/a = 1.25. 
The open circles and filled squares represent the results from the CPA and VCA, respectively. 
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Fig. 3. Partial Bloch spectral function of d x 2_ y 2 and d xy states at the T point in FeCo with 
x = 0.5 and c/a = 1.25: (a) results for the ordered structure (S = 1.0); (b) and (c) results for the 
disordered structure (S = 0.0) from the VCA and CPA, respectively. 



7/Q 



Appl. Phys. Express 



References 

1) T. Burkert, L. Nordstrom, O. Erikssonm, and O. Heinonen: Phys. Rev. Lett. 93 
(2004) 027203. 

T. Klemmer, D. Hoydick, H. Okumura, B. Zhang, and W. A. Sofia: Scripta Met. 
Mater. 33 (1995) 1792. 

G. Andersson, T. Burkert, P. Warnicke, M. Bjorck, B. Sanyal, C. Chacon, C. 
Zlotea, L. Nordstorm, P. Nordblad, and O. Eriksson: Phys. Rev. Lett. 96 (2006) 
037205. 

A. Winkelmann, M. Przybylski, F. Luo, Y. Shi, and J. Barthel: Phys. Rev. Lett. 
96 (2006) 257205. 

P. Warnicke, G. Andersson, M. Bjorck, J. Ferree, and P. Nordblad: J. Phys.: 
Condens. Matter 19 (2007) 226218. 

F. Luo, X. L. Fu, A. Winkelmann, and M. Przbylski: Appl. Phys. Lett. 91 (2007) 
262512. 

F. Yildiz, M. Przybylski, X.-D. Ma, and J. Kirschner: Phys. Rev. B 80 (2009) 
064415. 

C. Neise, S. Schonecker, M. Richterm K. Koepernik, and H. Eschrig: Phys. Status 
Solidi B 248 (2011) 2398. 

G. Hautier, A. Jain, and S. P. Ong: J. Mater. Sci. 47 (2012) 7317. 

B. Velicky, S. Kirkpatrick, and H. Ehrenreich: Phys. Rev. 175 (1968) 747. 
R. J. Elliott, J. A. Krumhansl, and P. L. Leath: Rev. Mod. Phys. 46 (1974) 465. 
J. B. Staunton, S. Ostanin, S. S. A. Razee, B. Gyorffy, L. Szunyogh, B. 
Ginatempo, and E. Bruno: J. Phys. C: Condens. Matter 16 (2004) S5623. 
Y. Kota and A. Sakuma: J. Appl. Phys. Ill (2012) 07A310. 
Y. Kota and A. Sakuma: J. Phys. Soc. Jpn. 81 (2012) 084705. 

H. L. Skriver: The LMTO Method (Springer, Berlin, 1984). 

I. Turek, V. Drchal, J. Kudrnovsky, M. Sob, and P. Weinberger: Electronic 
Structure of Disordered Alloys, Surface and Interfaces (Kluwer, Boston, 1997). 
I. Turek, V. Drchal, and J. Kudrnovsky: Philos. Mag. 88 (2008) 2787. 
R. M. Bozorth: Ferromagnetism (van Nostrand, New York, 1951). 
P. Bruno: Phys. Rev. B 39 (1989) 865. 

D. Wang, R. Wu, and A. J. Freeman: Phys. Rev. B 47 (1993) 14932. 
I. V. Solovyev, P. H. Dederichs, and I. Mertig: Phys. Rev. B 52 (1995) 13419. 



8/Q 



Appl. Phys. Express 



22) Y. Kota and A. Sakuma: in preparation for publication. 



Q/Q 



